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Abstract 

We calculate the rate of energy radiation by bremsstrahlung of neutrino pairs by 
baryons in neutron stars employing a transport model where neutrinos couple to 
baryons with spectral width. The coherence effects, which are included by computing 
the self energies with fully dressed propagators, lead to the Landau-Pomeranchuk- 
Migdal suppression of the neutrino emission rates in the soft neutrino and high 
temperature limit. A microscopic evaluation of the bremsstrahlung rate by neutrons 
is carried out using the Brueckner theory of nuclear matter at finite temperatures. 
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1 Introduction 



Neutrino pair bremsstrahlung by baryons, B\ + B 2 — > B\ + B 2 + v + 17, is 
among the dominant processes by which neutron stars radiate away their en- 
ergy during the first several thousand years of their evolution. The rates of 
these processes have been derived by Friman and Maxwell [1] and Voskresen- 
sky and Senatorov [2] employing the Fermi-liquid interaction combined with 
different prescriptions for pion exchange contribution to the nuclear matrix 
element. 



The spectator in the reaction above is needed to insure that the process is 
allowed kinematically.Q If the formation length of the neutrino radiation is of 

1 The neutrino pair bremsstrahlung via the reaction B — > B + v + V is inefficient as 
the energy conservation requires oj ~ qvp where q is the momentum transfer and vf 
is the baryon Fermi velocity, while for neutrinos in the time-like region oj > qc. In 
more formal terms the imaginary part of the quasi-particle single-loop polarization 
function vanishes in the low-temperature limit. 
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the order of the mean free path of a baryon, the role of the spectator is taken 
over by the medium because the baryon undergoes multiple scattering during 
the radiation. Then the reaction rate is subject to the Landau-Pomeranchuk- 
Migdal (LPM) suppression well-know from QED. At not too low temperatures 
the LPM suppression of the bremsstrahlung could be a factor in neutron stars, 
for the radiation frequency is of the order of temperature and could have a 
magnitude comparable to the baryon's collisional width. 

The LPM effect has been naturally recovered in quantum many-body theory 
in attempts to overcome the limitations of the quasi-particle approximation 
which led to infrared-divergent results in the soft limit. Knoll and Voskre- 
sensky, ref. [3], developed a general approach to the LPM effect using the 
Schwinger-Keldysh formulation of non-equilibrium field theory. The particle 
production rates in their approach are calculated from an expansion of the 
self-energies in terms of skeleton (or compact) diagrams with fully dressed 
particle propagators. Raffelt and Hannestad and Raffelt [4] have discussed 
in the supernova context the neutrino pair bremsstrahlung using the one-pion 
exchange interaction and including leading order corrections in particle width. 

In this Letter we shall adopt the formalism of ref. [3] to compute the rate of 
the neutrino pair bremsstrahlung from baryonic matter in neutron stars and 
supernovae, including the LPM effect. The spectral width of baryon propaga- 
tors, which enters the rates of the bremsstrahlung, is derived from the finite 
temperature Brueckner theory of nuclear matter. Friman and Maxwell's [1] 
result is adopted as the reference for comparison. 

The rates of neutrino pair bremsstrahlung in the Schwinger-Keldysh technique 
can be derived using the formalism of the optical theorem [5]. Here we follow 
a somewhat different path. It is known that the rates of the neutrino emission 
can be recovered from the Boltzmann equation for neutrinos [6]. For present 
purposes the Boltzmann equation, which deals only with on-shell particles, 
is not sufficient. We shall recover neutrino emission rates form the Kadanoff- 
Baym transport equation [7] which includes the collisional coupling of neutri- 
nos to baryons with a finite spectral width. As the dense cores of neutron star 
develop, apart from neutrons and protons, a large equilibrium fraction of hy- 
perons [10], our discussion shall refer generally to baryons. We treat neutrinos 
as relativistic massless Dirac fermions within the relativistic Kadanoff-Baym 
transport theory [8,9]. Our final expression relating the neutrino emissivity to 
the polarization of the medium agrees with the one derived by Voskresensky 
and Senatorov [5] from the optical theorem. 
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2 Formalism 



We start from the KadanofT-Baym transport equation for the neutrinos:^ 

{Y,9 Xfl S < (q,x)} + = a < (q,x)S > (q,x) - a> (g, x)S < (q, x), 

(1) 

where q = (g, go) and x are the neutrino four momentum and the center-of- 
mass space-time coordinate, respectively, 5 I> ' < and a > ' < (g, x) are the neutrino 
propagators and collisional self-energies, which are real-time contour ordered 
with fixed-time arguments on the lower (upper) and upper (lower) branches of 
the Schwinger-Keldysh contour, respectively. Treating neutrinos on-mass-shell, 
the propagators can be expressed in terms of the non-equilibrium distribution 
functions (Kadanoff-Baym ansatz [7]) 

5 (g - v v (q)) f v (q, x) - 6 (g + w v (q)) (1 - U(-q, x)) ] ; 

(2) 

the S > (q, x) propagator follows from eq. (2) via the interchange f u <-> (1 — /„). 
The collision integrals contain neutrino scattering, absorption and emission 
contributions of which only the latter one is relevant for the cold neutron 
stars. 

The (anti-)neutrino emissivity (the power of the energy radiated per volume 
unit) is obtained by multiplying the left-hand-side of the transport equation 
by the neutrino energy and integrating over an element of the four dimensional 
phase-space: 

= E ^ <,cm ® - J p* cm (^ "M, (3) 

where the reduced collision integral I^' em (q) (em stands for emission) origi- 
nates from energy integration over the range [0, oo] and corresponds to the 
neutrino branch of the spectrum, while Ip ,em (q) results from the integration 
over the range [— oo, 0] and corresponds to the anti- neutrino branch of the 
spectrum; the sum is over the neutrino flavors. The leading order Feynman 

2 Terms which do not contribute to the left-hand-side of the transport equation in 
the mean field approximation are dropped. 



S < (q,x) 



t . 1 I /"If 1 



3 



b 




Fig. 1. Dyson equation for (a) neutrinos (wavy lines) and (b) baryons (solid lines). 
Double and single lines correspond to dressed and bare propagators, respectively, 
the dashed lines correspond to the W, Z boson exchange. The light dots denote the 
neutrino weak interaction vertex. The dark dot stands for the particle-hole vertex 
renormalization in the baryon loop. The shaded block in the lower panel is the 
thermodynamic T-matrix in the lowest order Brueckner theory. 

diagrams in the expansion of the self-energies a >,< with respect to the weak 
neutrino-baryon coupling are shown in the Fig. la. The corresponding one- 
loop transport self-energies are read-off from the diagram 



-- > ' < ta.-) = /(gi/|^(^ 4 ta- 52 - 9 ) 

<r>:iS > ' < (q 2 ,x),rl x mX(-q,x), M 

where (q) are the Fourier transforms of the baryon polarization functions 
with time arguments fixed on the upper and lower branches of the contour. 
The baryon propagators in U^(q) are fully dressed; the Fig. lb displays the 
baryon propagator dressing due to the strong interaction within the lowest 
order Brueckner theory used in the numerical evaluation below. For small 
energy-momentum transfers involved in the problem, the weak interaction 
vertex, T^, can be replaced by the contact interaction: T M = ^G/2y / 2^ 7 M (1 — 
7 5 ), where G is the weak coupling constant. 

Inserting the self-energies and the propagators in the left-hand-side of the 
transport equation (1) for the neutrino branch, we find, e.g., for the first term 
which is the gain part of the collision integral 



(ft) 

5 (q 2 - UJ v (q 2 )) /„(g 2 , x) - 6 (q 02 + w„(ga)) (1 - /p(-?2, x)) 



4 



One may identify various contributions to the collision integral: the terms 
oc (1 — f v )f v and oc (1 — f v )(l — fp) in the gain integral correspond to the 
scattering-in and emission processes, respectively. Similarly the loss term, I>, 
which is obtained from eq. (5) via the simultaneous interchange f VjD <-> (1 — 
f v p) and replacement of IT> A by IT< A , contains the terms oc — /„) and 

oc f v f v corresponding to scattering-out and absorption, respectively. As well 
known, e.g. ref. [4], the neutrino mean free path exceeds the stellar radius 
and neutrinos leave the star without interactions except for several seconds 
after the neutron star birth. If neutrinos are untrapped, the only terms that 
contribute to the collision integral are the vv emission terms, one present in 
the neutrino gain integral (5) and the other in the anti-neutrino loss integral, 
Ip (which follows from eq. (5) by replacing II> A by n< A and changing the 
integration limits). 

Next we enforce the quasi-particle approximation in the collision integrals for 
neutrinos and anti-neutrinos by carrying out the q^ 2 and goi energy integra- 
tions. The gain and loss terms of the collision integrals are combined using 
the identities IT< A (g) = U> x (-q) = 2ig B (q )lmU^ x (q), where g B (q ) is the 
Bose function and U^ x (q) is the retarded component of the polarization func- 
tion. Substituting the electroweak vertex in eq. (3) and dropping the spatial 
argument, we find 



_ 2 ( G Vv r d3q2 f d3qi f d4q 

€vu \2V2J yJ {2^2u v {q 2 ) J (27r) 3 2^( ?1 ) J (2tt) 4 

(27r) 4 <5 3 (g 2 + qi - q)5(u u (q 2 ) + u v {qi) - go) 

M&) +u„(qi)\ g B (q )A^( qi ,q 2 )lmU^(q), (6) 

where A^g!,^) — Tr[7 M (l — 7 s ) ^i7 A (l — 7 5 ) jq 2 ] is the trace of the neu- 
trino currents. Here we made the approximation /„, C 1 appropriate for 
untrapped neutrinos. The final expression for the emissivity, eq. (6), recovered 
here from the transport equation, agrees with that derived from the non- 
equilibrium S- matrix expansion [5]. Analogous expressions for the scattering 
and absorption in equilibrium have been derived recently in refs. [11,12]. 

The Lorentz structure of the baryon polarization function H^q) is given by 
the trace of the baryon transition currents 

A m aG», q) = Tr[7 M (c y - ca75)(^ + m)~j X {cv - c A l^+ 4+m)\. 

The contraction of the neutrino trace with the non-relativistic reduction of 
X mA simplifies the factor k^ x (q u q 2 )lml\^ x (q) = %uj v (qi)uj v (q 2 )[c v L^(uj,q) + 
3c\L^(uj , q)] where the vector (Ly) and the axial vector (L^) one-loop polar- 
ization functions differ by the particle-hole vertex renormalization. The vector 
part corresponds to the scalar particle-hole vertex; the axial-vector part cor- 



5 



responds to the spin exchange interaction vertex. These renormalizations can 
be carried out [5] in the Fermi-liquid theory, using the quasi-particle approx- 
imation, and amount to replacing the weak-coupling parameters Cy and Ca 
with the effective ones (explicit expressions are given in ref. [5]). Eventually 
one is left with the factor 8c^(<fi)cov(<f 2 )(cy + 3c^)L^(cj,g) where 



L*(u>,$=2f- 



■ d 3 pde'de j g> (g e')g < (p + q, e) - g< (g e')g > (g + g, e) 
(2tt) 5 1 e' - e + w + iS 



(7) 



which is the driving term (and hence the subscript 0) in the resummation series 
for the particle-hole vertex; the parameters cy and ca here and below should 
be understood as the renormalized ones. The baryon propagators, g > ' < {p), are 
related to the spectral function a(p) and the Fermi distributions /jv(p) via the 
Kadanoff-Baym ansatz for baryons [7]: —ig < ip) = a(p) /jv(p), g < ip) —g > ij>) — 
ia(p), which is exact in the equilibrium limit. The spectral function, in turn, 
is given by 



a(u,p) = -21m [u - e p + ij(u;,p)/2] 1 , (8) 

where the quasi-particle energy e p = ep + ReS R (w,p)| a;=£p , e° p = p 2 /2m (i.e. the 
pole of the spectral function) and the inverse lifetime j(u,p) = — 2ImS i? (u;,p) 
(i.e. the width of the spectral function) are given in terms of the retarded 
baryon self-energy, T, R (p); our approximation for the latter quantity is depicted 
in Fig. lb. We note in passing that the imaginary part of the quasi-particle 
polarization function (the limit 'yip) — > 0) vanishes in the time-like region, 
9o > (fi m the low temperature limit. Eq. (7), therefore, is the leading order 
compact diagram contributing to the neutrino emission rate. 



3 Results 



Next we proceed to microscopic calculation of the self-energy and polarization 
function for the neutron matter, which is the dominant constituent of the 
dense stellar matter for densities n < 2n s (n s = 0.16 fm -3 ). The single particle 
energies and the width have been evaluated in the self-consistent lowest order 
Brueckner scheme at finite temperatures (ref. [13] and references therein). 
This scheme requires a self-consistent solution of the coupled equations for 
the thermodynamic (retarded) T^-matrix 
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Fig. 2. The real part of the on-shell self-energy and the half-width as a function of 
particle momentum at the saturation density n s for different temperatures; the zero 
temperature Fermi momentum is 1.7 fm . 

1 " fe>'' 2 ^(p,/)Gf(/,P,^)T^(/,p',P^),(9) 



+ 



7T 



E 



where a collectively denotes the quantum numbers (S, J, M) in a given par- 
tial wave, p and P are the magnitudes of the relative and total momentum 
respectively, V(p,p') is the bare nuclear interaction, is the angle averaged 
two-particle propagator 



Gf(p,P,w) 



l-f N (e(P/2+p)) l-f N (e{P/2-p)) 



4tt + -e(P/2-p) + i5 



(10) 
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with e(p) = Ep + ReT,(e p ,p) and the retarded self-energy 



= -E( 2J + !) dp'p'%?*{p,p'-p,p'-u + e{p'))f N {e{p')) ; 

la 

(11) 

here the lower case p is the magnitude of the proper momentum of a particle. 
The coupled equations (9) and (11) are subject to normalization to the total 
density at a given temperature. Fig. 2 displays the real part of the on-shell self 
energy and the half width of the spectral function as a function of the particle 
momentum for several values of the temperature at the saturation density n s . 
The numerical calculations were carried out using the Paris NN intercation 
keeping partial waves with L < 2. 

Let us next turn to the imaginary part of the polarization function, eq. (7). 
One of the integrations can be removed by the delta function; further only the 
angular integral can be performed analytically [3] . We substitute the resulting 
expression in eq. (6), carry out the neutrino phase-space integrals, and sum 
over their three flavors. The final result for the emissivity is 



?r 2 

" r 7 (4 + 34)m 2 T 7 J(T), (12) 



vv 5(2tt) 
where T is the temperature 



X(T)= dzz 5 g B (z) dy[f N (y)-f N (y + z)} dx 



\y-£(x)] +x(y,x) 2 /4 



arctan 



Z(x) + C(z)-y-z + 2y/xC(z) 
x(y + z,x)/2 

^(x) + ((z) -y-z- 2Jx((z) 



—arctan 



X(y + z,x)/2 



(13) 



with the quantities £(x) 
functions of x = £~/T, y 



~- e P /T, x(y,x) = j{uj,p)/T, ((y) = e q /T being 
lv/T, z = uj'/T; here e q is the recoil energy. Note 
that eq. (12) contains only universal couplings^ and constants as a prefactor, 



3 The values of parameters of weak interaction are cy = —1, ca = —(F + D) for 
neutrons, cy = 1 — 4sin 2 9yy, ca = F+D for protons, cy = — 2+4sin 2 0\y, ca = —F 
for S~ hyperons and zero for A hyperons; here F = 0.477±0.012, D = 0.756±0.011 
and sin 2 0\y = 0.23, where 6\y is the Weinberg angle. As noted above, the particle- 
hole vertex corrections renormalize these values of cy and ca- 
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all the nuclear many-body effects are absorbed in the integral X(T), which is 
valid at arbitrary densities and temperatures (in particular does not rely on 
any type of low-temperature expansion). At low-temperature our result can 
be simplified by expanding the spectral function (8) to the leading order in 7. 
Performing the same steps leading to eq. (12), we find in the first order 7 



(4 + 3c a) m*p F Z {p F )T 7 Jdz z 7 g B {z) x(z) 7- 



where p F is the baryon Fermi momentum, m* is the effective mass, and Z(p F ) 
is the wave-function renormalization; V stands for the principal value in- 
tegration. In this limit the width is essentially the reciprocal of the quasi- 
particle life-time in a Fermi-liquid, hence 7 cx T 2 [1 + (uj/2tcT) 2 ]. Substituting 
this in eq. (14) we recover the T 8 temperature dependence of the neutrino 
bremsstrahlung rate [1,2]. Note that according to ref. [1] the matrix element 
for the Fermi transition vanishes in the limit q —>■ in Born approximation. 
This cancellation occurs whenever the off-shell source particle propagation be- 
tween the strong and weak interactions is odd under time-reversal, which is 
true in the quasi-particle approximation neglecting the recoil term. We can 
not make an accessment on this cancellation, as the quantum interference 
diagrams are not included so far in our treatment. Should this cancellation 
indeed emerge as a consequence of the current conservation requirements/^] 
our numerical expression for the neutrino bremsstrahlung would overestimate 
the rate by a small factor Cy/3c\. 

A numerical estimate of the neutrino emissivity eq. (12) for the case of the 
neutron bremsstrahlung gives 



e v . = (7.56 x 10 18 erg cm" 3 s^ 1 ) T 9 7 J(T). (15) 

where Tg is the temperature in units 10 9 K. The triple integration in eq. (13) 
is carried out numerically after the functions £, x an d C are derived from the 
microscopic theory. Table 1 shows the results of the calculations of the integral 
X(T) and compares the neutrino emissivity, eq. (15), with the result of Friman 
and Maxwell[l]. 



4 Note that the arguments of ref. [1], perfectly valid in the limit discussed there, 
are not directly applicable to our discussion, as neither the particle width is an odd 
function of the frequency, nor the recoil term is neglected. 
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Table 1 



The upper number for a density-temperature pair is the value of the triple integral 
X, the lower number is the ratio of the eq. (15) to the result of Friman and Maxwell, 
eq. (66a). 



T (MeV) 


0.5 


2 


5 


10 


20 


30 


n s 


34.0 


19.3 


18.3 


8.6 


3.0 


0.4 




2.46 


0.35 


0.13 


3.1 x 10" 2 


5 x 10" 3 


5 x 10~ 4 


2n s 


33.0 


18.8 


17.64 


5.49 


3.84 


2.12 




2.38 


0.34 


0.127 


3.3 x 10~ 3 


7.0 x 10~ 3 


4.3 x 10~ 4 



4 Discussion 



In this work we recovered the relation [eq. (6)] between the neutrino emissivity 
and the polarization function of the nuclear medium from the Kadanoff-Baym 
transport theory. We derived the bremsstrahlung rates by computing the neu- 
trino transport self-energies which are driven by neutrino coupling to baryons 
with spectral width. The diagrammatic expansion for the polarization func- 
tion of the medium is constructed with fully dressed particle propagators, as 
originally suggested by Knoll and Voskresensky [3]. As a result the LPM ef- 
fect, i.e. the influence of the multiple scattering during the radiation process, is 
included in the rate of neutrino emission. We have kept only the leading order 
term in the expansion given by the one-loop polarization; more complicated 
diagrams are not considered at this stage. Already at the leading order the 
rates are well behaved in the soft neutrino limit and exhibit LPM suppression 
when 7/T > 1. Since this ratio at low temperatures is proportional to T the 
LPM suppression increases with increasing temperature; in our case the typi- 
cal scale for the onset of the LPM effect is T ~ 5 MeV. A comparison to the 
result of Friman and Maxwell [1] shows a strong suppression at relatively high 
temperatures and an enhancement by a factor of two in the low temperature 
limit .0 We believe that the latter low temperature enhancement is mainly due 
to the differences in the treatment of the nuclear interaction in the medium. 
For densities at and above the nuclear saturation density effects beyond those 
included in the Brueckner formalism, such as an enhancement of the strong 
interaction amplitude due to the in-medium modifications of the pion degrees 
of freedom, could compete with the LPM suppression, of. [2]. 



5 Note that the extrapolation of the Friman and Maxwell's result to high tempera- 
tures is limited to the assumption that the quasi-particle picture can still be applied 
and that the quasi-particles are lying close to their Fermi surface (the expansion 
parameter being the ratio of the temperature over the Fermi energy). 
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Another feature of the expansion above is that the quantum interference terms 
appear only in higher order diagrams. The latter correspond e. g. to particle- 
hole insertions which transfer momentum between the particle and the hole 
Green's functions [3]. It is well known that at T = inclusion of these inter- 
ference terms is required to fulfill the Ward-Takahashi identities [14]. At finite 
temperatures the interference corrections are reduced; an explicit evaluation 
at T 7^ has to be performed yet. 

This work has been supported by the Stichting voor Fundamenteel Onderzoek 
der Materie (FOM) with financial support from the Nederlandse Organisatie 
voor Wetenschappelijk Onderzoek (NWO). 
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